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Abstract 

The increasing interest in the Miiller density-matrix-functional theory has led us to a system- 
atic mathematical investigation of its properties. This functional is similar to the Hartree-Fock 
functional, but with a modified exchange term in which the square of the density matrix 7(x,x') 
is replaced by the square of 7"'^/^(x, x'). After an extensive introductory discussion of density- 
matrix- functional theory we show, among other things, that this functional is convex (unlike the 
HF functional) and that energy minimizing 7's have unique densities p(r), which is a physically 
desirable property often absent in HF theory. We show that minimizers exist if < Z, and derive 
various properties of the minimal energy and the corresponding minimizers. We also give a precise 
statement about the equation for the orbitals of 7, which is more complex than for HF theory. We 
state some open mathematical questions about the theory together with conjectured solutions. 

PACS numbers: 31.15.-p, 71.10.-w 
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I. INTRODUCTION 



The basic goal of density-functional theory is to express the energy of a quantum- 
mechanical state in terms only of its one-particle density p(r) and then to minimize the 
resulting functional (the 'density functional') with respect to p(r) (under the subsidiary 
condition that J^^ p{r)dr = N = number of electrons) in order to calculate the ground- 
state energy of the system, which could be an atom or a molecule or a solid. Although 
the first - and by far most used and important density functional in theory, computation, 
and mathematical investigation of multi-electron systems - is the Thomas-Fermi functional 
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(Lenz 2J]), strong interest in the subject was triggered by Hohenberg and Kohn [20]. We 
refer the readerjnterested in the recent developments to the books by Eschrig and Gross 



and Dreizler [17| and the review 



30|. 



While this program is possible in principal, experience has shown that it is far from easy 
to guess the appropriate functional - especially if one wants the functional to be universal 
and not simply 'tuned' to the particular kind of atom or molecule under investigation. There 
are also pitfalls connected with the admissible class of functions to use in the variational 
principle [2f 



Whereas the external potential energy can easily be expressed in terms of the one-particle 
density, it is not known how to express the kinetic energy and the interaction energy in terms 
of p(r). Going from density- to density-matrix-functional theory eliminates the first problem 
altogether, since all expectations of one-particle operators can be expressed in term of the 
one-particle density matrix. The density matrix analogue of the Hohenberg-Kohn density- 
functional program was established by Gilbert [l^. See also 25]. 

The most difficult component of the density-functional to estimate is the exchange- 
correlation energy (which we shall henceforth simply call exchange energy), and it is that 
energy that will concern us here. Owing to this and other difficulties, it has been the 
tendency recently to replace the energy as a functional of p(r) by a functional of the one- 
body density matrix, 7(x, x'). In this way it is hoped to have more flexibility and achieve, 
hopefully, more accurate answers. 

Fermions have spin and it is convenient to write a particle's coordinates as x = (r, a) for 
a pair consisting of a vector r in space and an integer a taking values from 1 to q. Here 
q is the number of spin states for the particles which - in the physical case of electrons - 
is equal to 2. (In nuclear physics one sometimes considers g = 4.) We shall, however, call 
the particles electrons. Similarly we write for any function / depending on space and spin 
variables ^ 

/ /(x)dx=^ / /(r,a)dr, (1) 

i.e., J dx indicates integration over the whole space and summation over all spin indices. 
This allows us to write the density matrix 7 as an operator on the Hilbert space of spinors 
ijj for which J ]'?/;(x)]^ c/x < 00. Its integral kernel is 7(x, x'). 
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The Schrodinger Hamiltonian we wish to consider is 

N 



i=l ^ ^ 



(2) 



where 

*■ z. 



w^E^ (3) 

is the Coulomb potential of K > 1 fixed nuclei acting on the N electrons. The j^^ nucleus 
has charge +Zje > and is located at some fixed point Rj G M^. We define the total nuclear 
charge hj Z = Ylf=i ^j- The electron-electron repulsion R is given by 

R= Ir-'-r.r^ (4) 

l<i<j<N 

If one is interested in minimizing over the nuclear positions Rj, one also has to take into 
account the nucleus-nucleus repulsion e^U, of course, which is given by 

U= J2 ZiZj\Ki-Rj\-\ (5) 

l<i<j<K 

Since we will not be concerned with this question but rather consider the nuclei to be fixed, 
we will not take this term into consideration here. 



A. Hartree-Fock Exchange Energy 

The best known density-matrix-functional associated with ([2]) is the Hartree-Fock func- 
tional 

^^"^(7) = ^ tr(-V^) - / K(r)p,(r)dr + e'D{p„ p,) - e'X{j) , (6) 
where ^^(r) = Yla=i 7(^' ^) ~ Yla=i ^' ^) particle density, 

0(P,,)-U I t^**'. (7) 

Jm? JR3 |r — r I 

and where the exchange term is (note the sign in (I6l)) 

1 f f |7(x,x'^i2 



X(7) = -y J ^^j-^d^d^'. (8) 

As is well known, this functional is the expectation value of if in a determinantal 
wavefunction made of orthonormal functions ipi 

v^(xi, X2, ... x;v) = (iV!)-V2det^.(x,)|5=i , (9) 
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in which case 

N 

7(x,x') = J]¥..(x)^,(x')* . (10) 

i=l 

It is also well known that any one-body density matrix 7 for fermions always has two prop- 
erties (in addition to the obvious requirement of self-adjointness, i.e., 7(x, x') = 7(x',x)*) 
which are necessary and sufficient to ensure that it comes from a normalized A^-body state 



satisfying the Pauli exclusion principle, see e.g.. 



3G, 



< 7 < 1 as an operator and tr7 = A^, (11) 

where tr denotes the trace = J (ix7(x, x) = sum of the eigenvalues of 7. A simple conse- 
quence of (fTTj) is that the spin-summed density matrix (tr^. 7) (r, r') = J2a 7(^' ^''^)' which 
acts on functions of space alone, satisfies 

< tro-7 < g as an operator and tr(tro-7) = A^. (12) 



The HF 7 in (jlOl) has A^ eigenvalues equal to 1, and the rest equal to 0, but one could 
ignore this feature and apply (El) to any 7 satisfying ( |TT1) . If we do this, then we can define 
the HF energy (for all A^ > 0) by 

E^^{N) = inf{fH^(7) : < 7 < 1, tr 7 = A^} . (13) 

7 

(We say 'infimum' in (fT3|) instead of 'minimum' because there may be no actual minimizer 
- as occurs when N ^ Z = Zj.) A HF energy minimizer does exist when A^ < Z + 1, 



at least, and possibly for larger A^'s as well |3J, |35 |. 

It is a fact 29|] (see also [2i]) that E^^{N) is the infimum over all 7's of the determinantal 
form fix 01) . i.e., the determinantal functions always win the competition in (IT^ . Therefore, 
E^^{N) > Eq{N), where Eq{N) is the true ground state energy of the Hamiltonian ([2]). 

Thus, the HF density-matrix-functional has the advantage of providing an upper bound 
to Eq, but it cannot do better than HF theory. We know, however, that this is often not 
very good, numerically, especially for dissociation energies. 

Another disadvantage of S^^ is that the energy minimizer 7^^ (if there is one) may 
not be unique although, in some cases, it is known to be unique (see ^| for the Dirac- 
Fock equations). In fact it follows from Hund's rule that in many cases the spatial part of 
the wave function has a non-zero angular momentum and cannot, therefore, be spherically 
symmetric. 



A third point to note is that in HF theory the electron Coulomb repulsion is modeled 
by D{p^, p-y) — X(7). This energy really should be Jj^g Jj^g |r — r'!^-*^ p*'^''(r, r')dYdv' , however, 
where p^^''(r,r') is the two-particle density, i.e., the spin summed diagonal part of the two- 
particle density matrix. In effect, one is replacing p*^^^(r,r') by G'*^^^(r, r') = ip^(r)p^(r') — 
\ X]cto-'=i |7(^>^')P- It is not possible for this G^^^ to be the two-body density of any state 
because that would require that Jj^g (^(^^(r, r')c?r' = ^^^pir). This condition fails unless 
the state is a HF state (because even the total integral is wrong, namely, jj G^'^^drdv' > 
N{N — l)/2 unless we have a HF state). 



B. Miiller's Square-Root Exchange- Correlation Energy 



There is an alternative to £^^^(7), which we will call S^{'~f) (Miiller 38|]). It replaces 
the operator 7 in X{'y) by 7^/^. This means the operator square root (note that 7 is self- 
adjoint and positive as an operator, so the square root is well defined). Thus, 7(x, x') = 
J (ix"7^/^(x, x")7^/^(x", x'). In terms of spectral representations, with eigenvalues Aj and 
orthonormal eigenfunctions (pi (the 'natural orbitals'). 



7(x,x') = 5^ A.y..(x)v,,(x')* and 7'/'(x,x') = ^ V,(x)^,(xO*. (14) 

1=1 i=l 

There is no simple formula for the calculation of 7^/^(x, x') in terms of 7(x, x'), unfortunately, 
but there is an integral representation, which we shall use later. 
Thus, 

^2 



2m 



tT{-V'-f) -e' [ V,{r)p,{r)dr + e^D{p„ p,) - e^X{^^''') , (15) 

JR3 



and 



E^{N) = inf{f'^^(7) : < 7 < 1, tr7 = A^} . 



(16) 



The functional £^^^(7) was introduced by Miiller 38| and was rederived by other methods 
3y_^Buijse and Baerends [5]. A similar functional was introduced by Goedecker and Umrigar 



ij], the chief difference being that [l4| attempts to remove an electron 'self-energy' by 



omitting certain diagonal terms that arise when ( |T6l) is explicitly written out using the 
expansion of 7 into its orbitals (fT4l) . In particular, quite analogous to density functional 
theory, explicit corrections terms have been added to correct the overestimate of binding 
energies using Miiller's functional (Gritsenko et al. 15|). 



From now on we will use atomic units, i.e., h = m = e = 1. To get some idea of the 



magnitudes involved we can look at hydrogen. Numerical computations [15|, Figure 6] and 
ISl . Figure 3.1] suggest that E^{1) ^ —0.525. This is to be compared with the true energy, 



-0.5. 

It might be wondered how Miiller's exchange energy compares to the old Dirac 
— / p-^(r)^/^(ir. As remarked after Lemma [21 and as found earlier by Cioslowski and Pernal 
j^, X(7^/^)can not be bounded by C / p^{r)'^^'^dr for any C. 



Miiller 38| also considered using 7P(x, x')7^^^(x', x) for some < p < 1 in place of 
|7"'^/^(x, x')p = 7"'^/^(x, x')7^/^(x', x), which satisfies the integral condition, but he decided 
to take p = 1/2 because this yields the smallest value of X, and hence the largest energy. 
(The proof is analogous to aPh^~^ + a^'^lf > for positive numbers a, b.) 

Miiller's functional f|T5|) has several advantages, the first of which is 
A.l. The quantity that effectively replaces p^'^\r,r') in the functional is now 

<? 



^p^(r)p^(r') - ^ J2 

a,a' = 

and this satisfies the correct integral condition 



|y/^(r,a, r',a')r, 



p,(r)p,(r')- Yl y/'(x,x')y/^(x',x) 

o-,cr' = l 



On the other hand, P7(r)P7(r') ~X]CTcr'=i 17^''^ (x, is not necessarily positive as a function 
of r, r', whereas the HF choice p^(r)p^(r') — X]cto-'=i |7(x, > (which is true for any 
positive semi-definite operator). This non-positivity is a source of some annoyance. In 
particular, it prevents the application of a standard method [s]]] for proving a bound on the 
maximum N. 

A. 2. A special choice of 7 is a HF type of 7, namely one in which all the Aj are or 
1. In this special case 7^'^^ = 7 and the value of the Miiller energy equals the HF energy. 
Thus, the Miiller functional is a generalization of the HF functional, and its energy satisfies 
E^{N) < E^^{N) (because, as we remarked above, the minimizers for the HF problem 
always have this projection property). 

Later, we shall propose that the quantity E^{N) = E^{N) + should be interpreted 
as the binding energy; it is not obvious that E^(N) satisfies such an inequality, however. 
Indeed, it does not, in general, as the hydrogen example shows (—0.525 + 1/8 > —0.5). 
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A. 3. The original Miiller functional seems to give good numerical results when few 
electrons are involved. Moreover, E^{N) appears to satisfy E^{N) < Eo{N) for all electron 
numbers A^, i.e., it is always a lower bound. We shall prove this inequality when = 2 in the 
last section. (Numerical accuracy of larger electron numbers seem to require appropriately 
modified functionals. We refer the reader interested on numerical results and improved 
density matrix functionals to the papers of Buijse and Baerends jsl, Staroverov and Scuseria 



45 1 . Herbert and Harriman 19|, Gritsenko et al. [15], Poater et al. [41j, Lathiotakis et al. 
22 1, and Helbig [3].) Since we are primarily interested in the structure of the underlying 
theory rather than numerical results, we concentrate on the unmodified original Miiller 
functional despite the above mentioned numerical deficiency for large electron number. The 
Miiller functional can be viewed as a prototype of density matrix functionals with simple 
structures, but which are potentially useable as the basis of more elaborate functionals, e.g., 



15|. 



C. Convexity and Some of its Uses 

A key observation about S^{'~f) is that it is a convex functional of 7. This means that 
for all < A < 1 and density matrices 71,72 (not necessarily with the same trace and not 
necessarily satisfying 7 < 1) 

£^(A7i + (1 - A)72) < Af ^(71) + (1 - A)^^(72) • (17) 

(Note that the convex combination A71 + (1 — A)72 satisfies the conditions in ( |T6l) if 71 and 
72 both satisfy the conditions.) The convexity is a bit surprising, given the minus sign in 
the exchange term of S^, and it will lead to several important theorems. One is that the 
electron density p^ylr) of the minimizer (if there is one) is the same for all minimizers with 
the same N, and hence that the density of an atom is always spherically symmetric. This 
contrasts sharply with HF theory, whose functional is not convex, and it can contradict 
the original Schrodinger theory (since an atom can have a nonzero angular momentum in 
its ground state). Also, the Dirac estimate for the exchange energy, — J p^l^ is not convex; 



it is concave, in 
Some writers 



'act! 

16=] regard the retention of symmetry as a desirable property for an approx- 



imate theory; one speaks of the "symmetry dilemma" of HF theory (which means that while 
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symmetry restriction of HF orbitals improves the overall symmetry it raises the minimum 
energy). Miiller theory has no symmetry dilemma! 

From another perspective the sphericity of an atom might be seen as a drawback since real 
atoms sometimes have a non-zero angular momentum, and such states are not spherically 
symmetric. Sphericity is not a drawback, in fact, since density-matrix-functional theory deals 
with density-matrices obtained from all A^— particle states, including mixed ones (because 
the only restriction we impose is that the eigenvalues of 7 lie between and 1, and this 
condition precisely defines the set of 7 obtained from the set of mixed states, not the set of 
pure states). In the case of atoms there is always a mixed state with spherical symmetry, 
namely the projection onto all the ground states, divided by the degeneracy. This is the 
state that one sees (in principle) when looking at an atom at zero temperature (Liiders' 



projection postulate |37|). 

A second consequence of convexity is that the energy E^{N) is always a convex function 



of A^, as it is in Thomas- Fermi theory, for example [29|, l36||. This means that as we add one 
electron at a time to our molecule, the (differential) binding energy steadily decreases. Such 
a property is not known to hold for the true Schrodinger energy Eo{N). 

The convexity of ^^^^(7) is not at all obvious. All the terms except — X(7^''^) are clearly 
convex. In fact, the term D{p^, p^) is strictly convex as a function of the density p^{r) (strict 
inequahty in f|T7|) when p^^ 7^ p^J since the Coulomb kernel |r — r'|~^ is positive definite. It 
is this strict convexity that implies the uniqueness of P7(r) when there is a minimizer. 

To show convexity of £^(7), therefore, we have to show concavity (like ( |T71) but with the 
inequality reversed) of the functional X(7^/^). First, we write |r — r'l^"*^ = J^Bx{r)*Bx(r')dX 
where A is in some parameter-space A. There are many ways to construct such a decompo- 



sition. One way is due to Fefferman and de la Llave pjj, which we shall use in the sequel, 
in which the functions Bx are all characteristic functions of balls in M.^ and A parametrizes 
their radii and centers. Another way is |r — r'|^^ = C J^s |r — z|^^|r' — z\^'^dz. Anyway, it 
suffices now to prove that J (ix(ix'7^/^(x, x')i?(r)*7^/^(x', x)i?(r') is concave in 7, for any 
fixed function B(r). We can write this in abstract operator form as tvy^^'^B'^'j^^'^B. The 
concavity of such functions of 7 was proved by Wigner and Yanase |47| in connection with 
a study of entropy. 

Convexity also holds for Miiller's general p functional, which we mentioned earlier. It 
uses 7^(x, x')7-'^^''(x', x) in the exchange term. The fact that tT'~fPB'''~f^~PB is concave for all 
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< J9 < 1 was proved in 28|] and plays a role in quantum information theory 39 1. 

Another important use of the convexity of ^^^(7) is to significantly simplify the question 
of the spin dependence of 7(r, cr, r', a'). For concreteness, let us assume the usual case of two 
spin states (g = 2), but the conclusion holds for any q. In the HF problem it is not obvious 
how 7 should depend on a, a' and usually one makes some standard a-priori assumption, 
such as that 7^^(r, cr, r', a') = 7|,-f(r, r', + 7l^l(r, r', )(5o-,j5o-/^j. In the Miiller case 

this problem does not arise. Note that the functional is invariant under simultaneous 
rotation of a and a' in spin-space. (This means that we regard 7 as a 2 x 2 matrix whose 
elements are function of r, r'. The spin rotation is then just a 2 x 2 unitary transformation of 
this matrix.) If we take any 7(r, a, r', a') and average it over all such simultaneous rotations 
we will obtain a new 7 whose energy S^{^) is at least as low as that of the original 7 (by 
convexity). But 7 is clearly spin-space rotation invariant, which means it must have the 
form 

7(r,a,r',a') = ^7(r,r')®I (18) 
where I is the 2x2 identity matrix. The subsidiary conditions become 

tr7 = y 7(r, r)dr = N and < 7 < 2 . (19) 

The change from 1 to 2 in fll9p is to be noted. Often 7 is called the spin-summed density 
matrix. 

The conclusion is that to get the correct minimum energy one can always restrict attention 
to the simpler, spin-independent 7, but with the revised conditions (fT9l) . This is a significant 
simplification relative to HF theory. In much of the sequel we utilize the formal notation x 
instead of r, but the reader should keep in mind that one can always assume that 7 has the 
form flTSl) and all spin summations become trivial. 

A question will arise: Although it is possible to choose 7 in the form f|T8|) . are there 
other possibilities? They will certainly exist if 7 is not unique, (but we conjecture that 
it is unique since its density is unique, as we said). Even if 7 is unique we still might 
have other possibilities, however, when N is small. For example, we could take 7(x, x') = 
7(r, r') X 5cr,|5o-',t' ^^is density matrix is bounded above by 1 only if 7 < 1 (not < 2). 
This situation can arise if is small, but we expect that it does not arise when A^ > 1. In 
any case, we show that, for large A^ and Z, 7 has at least one maximal eigenvalue, namely 
2 (see Prop. E]). 
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In short, it is likely that whatever the Miiller functional has to say about the energy, it 
probably has little to say, reliably, about the spin of the ground state. Unlike HF theory. 



we do not have to worry about spin here. JThis does not mean that 
better as concerns spin. Sometimes it is 



i'' theory is necessarily 
1 



and sometimes it is not 2 
In the atomic case £^^(7) is also rotationally invariant and we can apply the same logic 
used above for the spin to the simultaneous rotation of r, r' in 7(r,r'). The conclusion is 
that we may assume the following computationally useful representation: 

00 e 2/ -I- 1 °° 

7(r,r') = E E le{ry)YUOr)Yl^{e^,) = ^J2^,{r,r')P,{cose), (20) 

£=0 m=-e e=o 

where r = |r|,r' = |r'|. The Yi^rn are normalized spherical harmonics, 6r is the angle of 
the vector r, etc., Pe is the i^^ Legendre polynomial and 6 is the angle between r and r'. 
Another way to say this is that we can assume that the eigenfunctions of 7(r, r') are radial 
functions times spherical harmonics Ye^m and that the allowed m values occur with equal 
weight. This observation can simplify numerical computations. 

Any other symmetry can be treated in a similar way. For example, in the case of a solid 
there is translation invariance of the lattice of nuclei. By wrapping a large, finite piece of the 
lattice on a torus (periodic boundary conditions) we have a finite system with translation 
invariance and we can conclude, as above, that we can assume that 7(r, r') is also translation 
invariant, which means that 7(r, r'), viewed as a function of r + r' and r — r' is periodic in 
the variable r + r'. 

One obvious symmetry is complex conjugation {i ^ —i ) in the absence of a magnetic 
field. Convexity implies that in the spin- independent formulation any minimizing 7 must 
be real, as shown in Proposition [9] of Section IIVBI 



D. The Miiller Equations 

If the Miiller functional has a minimizing 7 (with tr 7 = A^) then this 7 satisfies an Euler 
equation. A minimizer does exist if A^ < Z as we show in Theorem [2l It is not altogether 
a trivial matter to write down an equation satisfied by a minimizing 7. Conversely, one can 
ask whether a 7 that satisfies this equation is necessarily a minimizer. We partly answer 
these questions in several ways. 
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1. Suppose that 7 satisfies tr7 = and that 7 minimizes £^^(7), i.e., £^^(7) = E^{N). 
Then we conclude (by definition of the minimum) that 

£^((l-t)7 + t7')>^''(7) (21) 

for all admissible 7' with tr7' = and for all < t < 1. Conversely, if tr7 = and if fl^Tj) 
is true for all such 7' and for some < t < 1 (with t possibly depending on 7') then 7 is a 
minimizer. Alternatively, it suffices to require that for all such 7' 

|£^((1 - t)7 + t7')l*=o = lim i [£^((1 - t)7 + W) - ^^(7)] > . (22) 

To see that 7 is a minimizer we exploit the convexity of the functional S^, which implies 
that ^^((1 - t)7 + t7') < (1 - t)£^{-f) + t£^{y), and hence, from or ([22]), that 
^^(7) < ^^(7')- (Note that the convexity also implies that S^{{1 — t)7 + ^7') is a convex 
function of t in the interval [0, 1], which, in turn, implies that the right derivative defined in 
( I22I) always exists.) 

To summarize, we say that the equation defining a minimizer is (l22l) (for all 7'). To make 
this more explicit we have to compute the derivative in fl22l) . 

2. The variational equations are most conveniently written down in terms of 7^^^, 
the square root of a minimizer. In Proposition \T0\ we will show that 7^/^(r, r') satisfies 
the following variational equation. Let ip^ denote the effective potential v^7(r) = Vc(r) — 
/ p../{r')\r — r'l^Mr', where p^{r) = X]o-7(^'^) ~ Xlo- / |7^^^(x, x')p(ix' denotes the particle 
density. Then 

f-i^l - iV?, - ^,(r) - ^,(r') - - 2i?\ 7^/^(x, x') = 2e.^.(x)^,(x')* (23) 

where fx < —1/8, Cj < and ipi{x) is an eigenfunction of 7^/^ with eigenvalue 1, i.e., 
J 7^/^(x, x')?/'j(x')(ix' = '?/'j(x) for all i. Note that the number of -ipiS corresponding to 
eigenvalue 1 is necessarily less than A^. 

Conversely, is it true that any 7^/^ satisfying < 7^^^^ < 1 (as an operator) and 
tr (7^''^)^ = tr7 = A^ which is a solution to ( |23l) under the constraints mentioned above, is 
a minimizer of £^^(7)? Unfortunately, we can answer this question affirmatively only if we 
know that the density p-y{r) does not vanish on a set of positive measure. Presumably such 
a vanishing does not occur, but we do not know how to prove this and leave it as an open 
problem. 
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3. As a practical matter it is the fact that 7 satisfies (123|) that is important because it 
gives us equations for the orbitals of 7. A minimizer 7 can be expanded in natural orbitals 

'ipji'x) as 

7(x,x')=^A,^,(x)^,-(x')* 

j 

with corresponding occupation numbers (eigenvalues) < Aj < 1. Then 7^/^(x, x') = 

^ 1/2 

\j ipj{'x)ipj{'x')* . Multiplying (!23l) by ipii'x') and integrating over x' yields an eigenvalue 
equation for the ipil'x), namely 

[(-iV^ - 7^/' + 7'/' (-iV^ - y.,)] 1^.) - [Z, + 2/xAp) |^,) = 2e.|^.) . (24) 

Here, is the operator with integral kernel 

Z,(x,x')=7'/'(x,xO|r-rr^ (25) 
Taking the product with {tpj \ , this implies, in particular, that 

{ipj\ - - >f-y\ipi) — j={ipj\Z^\ipi) = (/i + Ci) Sij. (26) 



(See also Pernal [40l] who derived - although merely on a formal level - similar equations 
for more general functionals). 

4. We shall show that 7 has no zero eigenvalues unless the density p-yir) vanishes 
identically on a set Q of positive measure. We do not expect such a set to exist but we 
do not know how to exclude this possibility. Any non-zero, square integrable function that 
vanishes identically outside is a zero eigenvalue eigenfunction of 7. In any case, there are 
no other zero eigenvalue eigenf unctions! 

Hence the orbitals ipji'x^) form a complete set in L^(]R^\i7). Formally, we can thus rewrite 
Eq. fl26|) as an eigenvalue equation for a linear operator on L^(]R^ \ Q). Let 

/7^ = -iV2-(/?^-X^, (27) 

where is the nonlocal exchange operator with matrix elements {ipilX^lipj) = (V^ + 

{ipi\Z^\ipj) . Alternatively, one can write 

1 f°° 1 1 
X^ = -/ Z^ y/^ds. (28) 

The variational equations are then 

H^\iPj) = fi\tp,) (29) 
14 



for all j with < Aj < 1, where /x < —1/8 is the chemical potential. Notice that all 
eigenvalues in fl29l) are identical, namely /i. 

In the subspace in which 7 has eigenvalue 1, which can only be finite dimensional since 
tr7 = A^, there is an orthonormal basis such that 



with all Cj < 0. The finite collection of numbers /i + ej constitutes all the eigenvalues of 
that are less than fi. 

The reason we say that fl2U]) and fl5U]) are formal is that the operator is only formally 
defined by (1271) . Both and are unbounded operators. Their sum is defined as a 
quadratic form (i.e., expectation values) but this form does not uniquely define the operator 
sum. If we knew that there are no zero eigenvalues then the set Q would be empty and 
would be defined as the usual Laplacian on M^, but if \ i7 has a boundary there 
are many extensions of with different boundary conditions, and this prevents the precise 
specification of (129]) and (1301) . There is no problem with the matrix elements in (1261) . however, 
since the ipi vanish on the boundary of \ Q. 

On the other hand (l23l) . which is an equation for the function 7^/^(x, x'), is true on the 
whole space. It is not necessary to impose any boundary conditions and is just the usual 
Laplacian - whether or not the set fl is empty. 

Surely fl is empty, in fact, and the practical quantum chemist can freely use (!29l) and 



E. Other Considerations about the Miiller Functional 

Let us conclude this introduction with a list of other significant questions about E^{N) 
and with statements about what we can prove rigorously. 

Ql. If there are no nuclei at all {K = 0), and if we try to minimize S^{'~f) (with tr7 = A^, 
however) it is clear that there will be no energy minimizing 7. There will, of course, be a 
minimizing sequence (i.e., a sequence 7„, n = 1, 2, .... such that £^^(7^) E^{N) as n — cxd. 
Such a sequence will tend to 'spread out' and get smaller and smaller as it spreads (always 
with tr7„ = A^). What, then, is E^{N)7 We prove that it is exactly given by 



(30) 



N/8 



when all Zj = 0. 



(31) 
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(If the units are included the energy is —{me'^/8h'^)N.) A similar calculation in the context 



of the homogeneous electron gas was done by Cioslowski and Pernal [6|] 

This situation is reminiscent of Thomas- Fermi-Dirac theory 29|] where, in the absence of 
nuclei, the energy equals —{const.)N. This negative energy comes from balancing the kinetic 
energy against the negative exchange. In such a case it is convenient to add + (const.) tr 7 to 
£^^(7) (with (const.) = 1/8 in our case) in order that E^{N) = when there are no nuclei. 

Another way to say this is that the energy, —1/8, is the self-energy of a particle in this 
theory. It has no physical or chemical meaning but we have to pay attention to it. It is the 
quantity 

E^{N) = E^{N) + - (32) 

8 

that might properly be regarded as the energy of electrons in the presence of the nuclei, 
i.e., —E^(N) is the physical binding (or dissociation) energy. We do not insist on this 
interpretation, however. On the other hand, if we are interested in the binding energy with 
fixed (e.g., the binding energy of two atoms to form a molecule) then it makes no difference 
whether we use the difference of E^{N) or E^{N). 

The motivation here is to ensure that the ground state energy of free electrons is zero. 
This can be compared with the formulation in [l^ in which the 'self-energy' correction is 
obtained by omitting certain diagonal terms in the energy (when the energy is written in 
terms of the orbitals of 7). This procedure does not have a natural physical interpretation 
and, more importantly, does not appear to give the zero energy condition for free electrons. 

This consideration leads us to the functional 

^(7)=^^(7) + ^tr7 (33) 

and its corresponding infimum E^(N). Note that ^(7) is also a convex functional of 7 
since the new term tr7/8 is linear, and hence convex. Likewise, E^{N) is a convex function 
of N. 

Having added this term, and with nuclei present, E^(N) will qualitatively look like the 
Thomas- Fermi energy, ^^^(A^). That is, E^{0) = and E^{N) decreases monotonically, 
and with non-decreasing derivative, as A^ increases 36|], 29|, Fig. 1]. It is bounded below, 
that is, 

E^{N) > E^ioo), (34) 
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where E^{oo) is some finite, negative constant. We shall prove this here. These features 
are displayed schematically in Fig. 1. 

There is another feature of E'^^{N) that we believe to be true for E^{N), but leave 
as an open question. At a certain critical value, Nc, of the electron number E^^{N) stops 
decreasing and becomes constant for all > A^^^. When N > the excess charge N — Nc 
just leaks off to infinity. In TF theory A^ is the neutrality point Z = ^ Zj, but this need 
not be so in other theories. In the original Schrodinger theory ([2]) A^ is greater than Z for 
many atoms (since stable, negative ions exist) but we know it is less than 2Z + 1 3l|]. In 



the Thomas- Fermi- Weizsacker theory, Nc is approximately Z + {const. 
know how to prove that there is a finite Nc for E^{N), but we believe there is one. 



29j. We do not 




E(N) 



FIG. 1: Schematic diagram of the energy dependence on the particle number A^. The lower, 
dashed curve is the Miiller energy E^{N) and the upper, solid curve is E^(N) = E^(N) + 
A^/8, in which the 'self-energy' —N/8 has been subtracted. Beyond the value Nc each curve 
is linear, whereas for N < Nc each is strictly convex and there is an energy minimizing 
density matrix. 

Q2. The main problem that has to be addressed is whether or not there is a 7 that 
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minimizes ^(7) in f|T6l) . If A^c < 00 we know that there is no minimizer when N > N^, so 
we obviously do not expect to prove the existence of a minimizer for all A^. 

The way around this problem, as used in [36|], for example, is to consider the relaxed 
problem 

E^{N) = inf{^(7) : < 7 < 1, tr 7 < A^} . (35) 

7 

The relaxation of the number condition allows electrons to move to infinity in case is 
larger than the maximal number of electrons that can be bound. In Proposition H] we show 
that E^{N) = E^{N) for all A^. 

The difference is that while the problem may not have a minimizer we prove that the 



E^ problem (l35ll has a minimizer for a. 
than the analogous proof in TF theory 



1 A^. The proof is more complicated in several ways 
29 . 3^. A minimizer, which we can call 7<(A^), will 



have some particle number tr7<(A^) = A^< < A^. It then follows from standard arguments 
using convexity (and strict convexity of D{p,p)) that the following is true, as displayed in 
Fig. 1 : 

If A^< < A^ then 7<(A^) = 7<(A^<) and E^{N) = (constant) = E^{N<), i.e., the original 
problem ( fT6|) has no minimizer. 

If A^< = A^ then 7<(A^) is also a minimizer for the original problem f|T6|) . That is, the 
relaxed problem and the original problem give the same minimizer and the same energy. In 
this case, E^{N) < E^{N') for all A^' < A^. The largest A^ with this property is equal to 

It might occur to the reader that nothing said so far precludes the possibility that Nc = 0, 
but this is not so. We prove that Nc> Z = total nuclear charge. 

Q3. How many orbitals are contained in a minimizing 7? We shall prove that 7 has 
infinitely many positive eigenvalues. This feature also holds for the full Schrodinger theory 
(Friesecke 12 1 and Lewin ^^), whereas there are only A^ in HF theory. We believe that 
7 has no zero eigenvalues (in the 'spin-summed' version), but cannot prove this. In other 
words, we believe that the eigenfunctions belonging to the nonzero eigenvalues span Hilbert 
space (they form a complete set). We can, however, prove that the eigenfunctions of the 
spin-summed 7 are a complete set on the support of p^{r), namely on the set of r G M'^ for 
which p-y{r) > 0. Presumably, this is the whole of M^. 

This introduction is long, but we hope it serves to clarify our goals and results, since the 
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rest of the paper is unavoidably technical. 
F. Open Problems 

For the reader's convenience we give a brief summary of some of the open problems raised 
by this work, some of which are discussed at various places in this paper. 

1. What is the critical value of the total electron charge, N^, beyond which there is no 
energy minimizing 7 and the energy E^(N) is constant? Is N^. finite and can one 
give upper and lower bounds to it? In particular is A^^^ > Z, i.e., can negative ions 
exist? (We prove Nc > Z and we prove that E^{N) is bounded below, for all A^, by 
a Z-dependent constant.) 

2. Is E^{N) < the true Schrodinger ground state energy? (We prove this for N = 2.) 
Can anything be said, in this regard, about E^{N) = E^{N) + N/87 

3. Is the spin-summed energy minimizing 7 unique? (We prove that all minimizers have 
the same density p(r), however.) 

4. Is the domain on which the unique p(r) > equal to the whole of (except, possibly, 
for sets of measure zero)? If so, this would imply that the spin-summed 7 does not 
have a zero eigenvalue. 

5. What are the qualitative properties of the density p(r)? How does it fall off for large 
|r|? What is its behavior near the nuclei? 

6. In this theory do atoms bind to form molecules? (Recall that there is no binding in 
Thomas- Fermi theory {sgI.) 

II. THE CASE Z = 

As noted in the introduction the energy of free electrons E^{N) is not zero but is pro- 
portional to A^. To be precise, E^{N) = —N/8 (in atomic units) when there are no nuclei, 
and comes about from the negative exchange energy — X(7^/^). This negative energy could 
be —00 were it not for the positive kinetic energy, which controls it and leads to a finite 
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result. We shall prove that the direct Coulomb repulsion term, D{pj,pj) plays no role here 
because it is quadratic in 7, whereas the terms we are concerned with are homogeneous of 
order 1. We would get —N/8 even if we omitted the direct term. Similarly, the value —N/8 
is independent of the number of spin states q. Moreover, the assumption 7 < 1 is not needed 
in the proof. 

In this section, Z = Zj = 0, and we are considering the functional 

£^(7) ^ tr(-iV%) + Dip,, p,) - X(7^/^) (36) 

and the minimal energy E^{N) in (|T6l) . We also consider the relaxed energy E^[N) for 
which, in analogy with (l35l) . the condition tr7 = is replaced by tr7 < A^. 

We always assume that (— + 1)^/^7^/^ g (5^, the set of Hilbert-Schmidt operators, so 
tr((l-V^)7) = / / c?xc?x' (|V7^/^(x,x')P + |7^/^(x, x')P) < 00. We use the usual notation 
for L^-norms, namely 

i/p 



p 



J |/(x)rrfx^ ' and ||/|U = sup{|/(x)|}. 
Proposition 1. If Z = 0, then for any N > 0, 

E^{N) = E^{N) = -N/8 (37) 

and there is no minimizing 7. 

Proof. Lower bound : We use the lower semi-boundedness of the hydrogenic Hamiltonian 
(i.e., for an imaginary nucleus with Z = 1/2, located at r') 

-iV?-(2|r-r'|)-i>-| (38) 

for all r' G M^, together with the fact that D{p,,p,) > to get 

^^(7) >lll (|Vry/^(x,x')r - '^'|^^"'"|^'' ) Cixdx' 

>-- [[ |7i/2(x,x')|'rfxc/x' = -itr7. 



This proves the lower bound on E^{N) and E^{N). 

To prove the non-existence of a minimizer we denote by g{r — r') the ground state of 
-Vr - |r - r'l"^ i.e., 

g{r-r') = n-'/^e-^'-'\ (39) 
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and note that the inequahty < in (138|) is strict (i.e., it is >), except for multiples of the 
function g(r — r'). Hence the above lower bound on £^{'j) is strict unless 7^/^(x, x') = 
Co-a-'(r')5f(r — r'). By self-adjointness, Cg-o-' has to be a constant, and since 7 G the set of 
trace class operators, Cg-o-' = 0. But this means that there exists no minimizer. 
Upper bound: We define a trial density matrix 7 by defining its square root: 

7^/^(x, x') = xirYgir - r')x(r') q-'/X,.'- (40) 

Here, g is the same as in (!39l) and x is a smooth function which will be specified later. Note 
that this definition makes sense, since the operator whose kernel is given on the right side of 
(HOl) is non-negative. This follows from the positivity of ^, the Fourier transform of g, given 

by 



^(P) 



23/2 



TT (l+|p 

An easy calculation shows that 



2\2' 



tr(-Vh) = {\xir)W{r')\'i~Vlg{r - T'))gir - r') + \Vx{r)\'g{r - r')^ \x{r')\') dv dv' . 
Using the eigenvalue equation for g one finds 

tr(-Vr')7 = 2X(7^/^)-itr7+ // | Vx(r) l^r - r')^ |x(r') T (ir ^ir' . 



4 

The upper bound will follow from this if we can find functions xl (where L is some free 
parameter) such that for 72, defined via xl, 

7l < 1, as an operator, tr72, N, (41) 

\VxL{r)\'g{r-r')' Ixdr^ dr dr' ^ 0, and D{p,^,p,J ^ (42) 

as L — s> 00. We shall choose xl of the form XLi^) = '^^^x(r/i^) for a fixed smooth function 
X > satisfying \\x\\t = N. 

We note that for any function ip (and with — denoting the Fourier transform) 



(^,7f V^) = i2nf'J ?(p)|(SO(p)r^P< (2vr)^/1?||oo||XL||Lll^ll2, 



which is less than or equal to \\ip\\2 for L large, since ||xl||oo 0. This implies the first 
condition in (j4T|) . To check the second one, we write 



tr7L = (2vr)3/^y" (?) (p) | ('^i) (p) ^ dp. 
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Now |(xi)(p)P = -^^l(x^)(-^p)P) which converges to N6{p) as L ^ oo (recall that ||xil4 = 
N). Therefore 

tr7L ^ (27r)3/2(^(0)A^ = A^. 
To check conditions (H2|) we estimate (again using that \\g\\2 = 1), 

II \VxLir)\'gir-rrXL{rrdrdr'<\\xL\\l / l^xdr)]' dr = L-'llxWloW^xf ■ 
Moreover, 

where v^L(r) = J g^{L{r — r'))x^{r') dr' . Since v^L(r) — s> x^(r) as L ^ oo, we conclude 
that D{p^j^, pyj^) = L^^D[x*] + o(L^^) by dominated convergence. 

Hence ( H2i) holds, and the proof is complete. ■ 

Remark: One might ask whether X(7^/^) can be bounded from above in terms of the 
usual Dirac type estimate for the exchange energy, J py{r)^^^ dr (cf. 29(]). However, this 
is not the the following example shows: define 7l, as in the proof Proposition [H 

by 7]^^^(x, x') = L^'^/^x(r/L)5f(r — r')x(r'/L)g^^/^5o-_o-', and carry out calculations similar to 
those done above. We find that 

|2 

2|r|' 



i/2n II ||4 / 9{r)t , 



Pyd^y^^dr^ L ^WxWfJ/l , 
P7lW dx 11x11 4 • 

Hence a bound in terms of the 4/3-norm can not hold. This example can be traced back to 

n 

Cioslowski and Pernal [6|. 



III. MINIMIZER IN THE CASE Z >0 



We return here, and in the remainder of this paper, to the general case in which all 
Zj > 0. We investigate the functional £^ in (155]) and the corresponding relaxed minimization 
problem given in ([35]) • Our goal is to show that there is an energy minimizing 7 for this 
problem and that its trace is tr7 = whenever N < Z = Zj. The main result of this 
section is contained in the following two theorems, whose elaborate proof will be given in 
several parts. 
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Theorem 1. For any Z > and N > one has E^{N) < and the infimum fl35l) is 
attained. 

As explained in the introduction, we do not know how to prove that the minimizer is 
unique. The strict convexity of the direct energy D{p^,p^), however, does imply that all 
minimizing 7's have the same (spin summed) density p-yir). 

Theorem 2. Assume that N < Z . Then a minimizer of (135|) has trace N. 

In particular, this result implies that in the original problem (fT6l) the infimum is achieved 
in case N < Z. The critical number Nc mentioned in the introduction is thus at least Z. 



A. Proof of Theorem [T] 

By Proposition [H the functional ^(7) is non- negative, ii Z = 0. By using a trial density 
matrix, we will first show that it assumes negative values as soon as Z is positive. 

Lemma 1. For any Z > and N > one has E^{N) < 0. 

Proof. Without loss of generality we may assume that there is only one nucleus of charge Z 
located at the origin r = 0. We use the same family 7^ of trial density matrices as in the 
proof of the upper bound in Proposition [1] Using the same estimates, we have 

S^{^l) = -Ztr |r|-SL + Y^lx'] + o{L-') as L ^ 00. (43) 

Since J 5f^(L(r — r'))x^(r') dr' X^(r), we have tr |r|^-'^7L = L^^ J |r|^-'^x^(r) dr + o{L^^) . 
Hence, 

£^{1l) = (^-Z J |rrY(r)rfr + D[x1^ +o(L-i) as L ^ 00. (44) 

For Z > and = ||xll4 small enough, the first term in brackets can clearly be made 
negative by an appropriate choice of x- This shows that E^{N) < for small N, and hence 
for all N. m 

Proposition 2. Let Z > and N > 0. There exists a minimizing sequence 7^ for ( l35il which 
converges in the space of trace-class operators, i.e., there is a'j such that ti |7j — 7I 0. 

Before giving the proof of this proposition, we collect some useful auxiliary material. 
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Lemma 2. For every e > 

'{|r-r'|<e} 

and 



[ [ b'^H^,^')? ^ 4^ tr(-V')7 (45) 



X(7^/')<|tr(-V^)7 + ^tr7. (46) 
Proof. The first inequality can be easily deduced from Hardy's inequality, which states that 

For the second inequality, we use the well known expression for the ground state energy 
of the hydrogen atom, namely. 



from which it follows (with z = 2/e) that for every x' 

1 | |7nx,xOp ^^^.y- |vy/2(^,^')|2rf^+^ I |y/2(^,^')|2rf^. (49) 

The lemma follows by integrating over x'. ■ 
Lemma 3. Let x(r) satisfy |x(r)| < 1- Then 

Proof. For convenience we introduce the characteristic function of a ball of radius r centered 
at z 

{1 |r — z| < r 
(50) 
|r — z| > r. 

Writing the Coulomb kernel as 

I r r ^ , . „ , . dr 



5.,.(r)5,,,(r')rfz— (51) 



vr ./o ./ro3 



(Fefferman and de la Llave [ll|), we get 

1 /"°° /" dr 
Xi^) = 7r tri6B,,rSB,,r) dz—. (52) 

It follows from |x| < 1 and the monotonicity of the operator square root that 
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Hence 

tr {ri"\B,,rX*l"\B,,r) < tr {{rixf'^B,,r{x*l"\f'^B,,) . 
The assertion follows now from fl52|) . ■ 

Proof of Proposition We choose an arbitrary minimizing sequence 7j for (!35l) and, after 
passing to a subsequence (if necessary), assume that tr7j —>■ N & [0,A^]. It follows from 
W> and the hydrogen bound, tr Zk\r - Rfcl'^ < (Zke/AZ) tr(-V2)7 + (ZkZ/e) tr 7 that 

^(1 - e) tr(-V^)7, < ^(7,) + V' + 1/4) tr 7,. (53) 



Hence the sequence 



V2 + 1)1/2 ^, (^_y2 _^ iy/2 bounded in and, by the Banach- 



Alaoglu theorem (see [33|]) there exists a 7 such that, after passing to a subsequence (if 
necessary), tr K'jj — > tr K'j for any operator K such that (— + l)^^^'^K{—V'^ + 1)~^/^ 
is compact. This compactness condition is satisfied if K is simply multiplication by some 
function / G Lp(M'^) for some 3/2 < p < 00 (see 42, section 13.4]). In this case we have 
that 

j /(r)p,^.(r)dr = tr/7,-^tr/7 = j /(r)p,(r)dr. (54) 

In particular, we can take / in ( 15^ to be the Coulomb potential since this potential can be 
written as the sum of two functions, one of which is in Lp(R3) and the other in L'^{R^) with 
3/2 < p < 3 and 3 < g < 00. 

Note that < 7 < 1 and, by the lower semicontinuity of the 6^-norm, 

M = tr 7 < liminf tr'jj = N < N. 

We claim that 7^0 (and hence M > 0). Indeed, by Proposition [1] one has ^(7^) < —e 
for some e > and all sufficiently large j. Hence tr V"c7j > e and by ( !54l) also trV^7 > e. 

Clearly, 7^- ^ 7 in the sense of weak operator convergence. If M = A^, then also tr7j — »• 
tr7, and thus 7^ — 7 in ©""^ (see Theorem A. 6 in \4^) and we are done. 

We are thus left with the case M < N. Our strategy will be to construct a minimizing 
sequence 7° out of the 7^- which converges to 7 in We choose a quadratic partition 
of unity, (x°)^ + {x^Y = 1? where x° is a smooth, symmetric decreasing function with 
;^0(0) = 1, x°(r) < 1 if |r| > and x°(r) = if |r| > 2. For fixed j, tr{x°{r/R)f-fj is a 
continuous function of R which increases from to tr7j. If we restrict ourselves to large 
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j, then tr7j > M and we can choose an Rj such that tr(x''(r/i?j))^7j = M. We write 
X^^r) = X'ir/Rj) and = X,"7,X," for = 0, 1. 

We claim i?j — 00. To see this, assume the contrary, namely that there is a subsequence 
that converges to some R < 00. Then, for this subsequence, X°(r)^ —>■ x*^(r/i?)^ strongly in 
any L^. Since p^. weakly in for 1 < p < 3 by fl54p . one has 



But, by definition, the left side is independent of j and equals / 

J p^{r) dr. This is a contradiction, since < 1 almost everywhere and 7^0. 

Therefore limj^oo Rj = 00 ■ We note that 7° ^ 7 in the sense of weak operator conver- 
gence. (It suffices to check the weak convergence on functions of compact support, since the 

remain uniformly bounded.) By construction, tr7° = tr7, so that 7° — 7 in &^ (again 
by Theorem A. 6 in 4^) and it remains to prove that 7^ is a minimizing sequence. 



For the kinetic energy we use the IMS formula 

tr(-VS) = tr(-V\°) +tr(-VSj) -tr[(|Vx,T + IVx^Ht,]- 
Since Rj ^ 00, one has |||VXjP + |Vx]Pl|oo and therefore 

tr(-V')7, = tr(-V')7° + tr(-V')7) + o(l). (55) 
For the attraction term we use again that Rj 00, so tr |r — Rfc|^-'^7j and 

tr |r - Rk\S = tr |r - Kk\~^^^ + o{l). (56) 
For the repulsion term we use that p^o < p-y. pointwise and get 

D (p^^. ,p^^)>D (p^o , p^o ) . (57) 
Finally, we turn to the exchange term, which we write as 

Xilf) = Xi^jfx^) + X{x]lTx]) + 2X(x?7;/^xj). 
We shall show that 

X{lf) < X((7°)^/^) + X{{^]y/') + o(l). (58) 
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It follows from Lemma [3] that X{Xj'yy'^Xj) ^ ^((7j)^^^)- To show that the off-diagonal 
term tends to zero we decompose, for any e > 0, 

J J{\r-r'\<e/2} A\r - V \ 

+ 11 — —. -p-^ axdx. 

ij{|r-r'|>£/2} 2|r — r I 

The term with the singularity is controlled by (HSl) . 

ff |x°(r)7/^(x,xOx](rOP ^ ^ ' ^ W Y72^ o 1/2, 1,2 1/2 

^^{|r-r'|<£/2} /:|r-r| 

<5tr(-V2)x;7,-X°- 

This can be made arbitrarily small be choosing e small. We pick some 5 > and decompose 
the term without singularity into two pieces, depending on whether |r'| < 5Rj or not. In 
the first case we estimate^ 

// lx;(r)7f(x.xOx;C-)P^^^^, ^^^j 

J J{\r-r'\>e/2,\r'\<5Rj} 2|r - r'| 

<e'' ff |7]/'(x,x')x](r')Nxrfx' 

J J{\r'\<SRj} 
= e-Hi X{\r\<SR,}{x]Ylj 
< ^''^N\\x{\r\<5R,}X]\\lo- 

Since x^ is smooth with x^(0) = 0, the supremum-norm of the function X{\r\<SRj}X] (which 
is independent of Rj by scaling) can be made arbitrarily small by choosing 6 small. Hence 
the double integral ( l59l) can be made arbitrarily small. 

In the complementary region one may argue as follows. We pick some A and choose j so 

[1] The following two paragraphs slightly differ from the published version in Phys. Rev. A 76 (2007), 
052517. We are grateful to M. Tiefenbeck for pointing out an error at this point of the proof. 
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large that Rj > S ^A. By estimating |r — r'| > SRj — A if \r\ < A and |r'| > SRj, we obtain 



|x°(r)7j/'(x,x')x](r')P 



dyid'x' 



{|r-r'|>e/2, |r'|><5_Rj} 2|r — r'| 



{|r-r'|>e/2,|r|>A} 2|r — r'| J y{|r|<A, jr' |>5_R,j } 2|r — r'| 

<e-^ II X?(r)'|7]^'(x,x')prfxrfx' + (2(5i?,- - / / |7j/'(x, x')P rfxrfx' 

JJ{|r|>A} 



-Urx{|r|>A}7,° + (2(5i?, - A))-4r7,-. 



Since 7° — > 7 in G^, one has tr X{|r|>yi}7j ~^ trX{|r|>A}7- This can be made arbitrarily small 
by choosing A large. Since Rj 00, the term {2{5Rj - A))-^ii 7j converges to 0. This 
proves fl58|) . 

Collecting ([SSD-(ISHD we find that 

^(7,) > ^(7°) + (-^trV^7j -X(7j) + ^tr7l) +o(l). 

We have shown in the proof of Proposition [1] that the term in brackets is non-negative. 
Hence 

liminf ^(7,) > liminf ^(7°), 
which shows that 7° is a minimizing sequence. This concludes the proof. ■ 
Proposition 3. Let 7^ 7 m (5"'^. Then 

liminf^(7,) >^(7). (60) 

Proof. The bound fl53|) shows that E = liminfj_^oo ^^{ij) > ~oo. Moreover, we may assume 
that < 00, for otherwise there is nothing to prove. After passing to a subsequence (if 
necessary), we may assume that S^{"yj) E. As in the proof of Proposition [2] there exists 
a 7 such that, after passing to a subsequence if necessary, ti K'jj — > tr_ft'7 for any operator 
K such that (— + l)~^/^-ft'(— + 1)~^/^ is compact. In particular, (l5il) holds. By weak 
lower-semicontinuity we infer that 

tr (-iV^ + 1/8) 7 < liminf tr (-^V^ + 1/8) 7,-. (61) 

Now we turn to the repulsion term. Since D{p^.,p^.) is bounded we may, passing to 
a subsequence (if necessary), assume that p^^ converges weakly to some p with respect to 

28 



the D-scalar product. With the help of (15^ one concludes that p = p-y. Weak lower- 
semicontinuity with respect to the D-norm implies that 



D{p^,p^) < lim inf D{p^ , p^. 



(62) 



The continuity of the attraction term follows from ( l54l) . since |r|^^ G L^/^ + for p > 3, 
therefore 

lim tr Vc'-yj = tr V^7 . (63) 

Finally, we prove continuity of the exchange term. Similarly as in the proof of Proposition 
[2] we decompose, for any e > 0, 

|7;/'(x,x')P + |7^/'(x,x'^|2 



{lr-r'|<£/2} 
+ 



2r-r' 



(ix dx' 



7^ (x,x')p-|7^/'(x,x' 



2|r-r'| 



(ix (ix' 



'{|r-r'|>£/2} 

According to Lemma [2] the term involving the singularity is bounded by £:tr(— V^)(7j + 7), 
which can be made arbitrarily small (recall that tr[— V^(7j + 7)] is bounded). To treat the 
term without the singularity we use the fact that the mapping K h-^ ji^'l^/^ is continuous 



from to & (see Example 2 after Theorem 2.21 in 4^). Hence 7^^^ 7^/^ in and 
we can bound 



< 



{|r-r'|>e/2} 



1/2 



ii7r(x,x') 


2 _ 


y/2(x,x')n 


2 


r — 


r' 





(ix (ix' 



1/2/ 



7/ ix,xj - 7 ' ^x,x 



c/xrfx' / ' " ^ 

'{|r-r'|>e/2} 



(ix' (ix' 



(| 7f (x,x')| + |7^/^(x,x')|)^ 
4 1 r — r' I ^ 

<||7f -y/^ll^2.-^tr(7,+7). 



The first factor tends to zero by the convergence of 7)^^ mentioned before, and the second 



one remains bounded. Hence we have proved that 

lim Xi^y^) = X(7i/2). 



(64) 



By collecting (!6Tl) - (l64l) we arrive at (160 



Proof of Theorem [II According to Proposition [21 there exists a minimizing sequence that 
converges strongly to some 7. By Proposition [31 this 7 is a minimizer of ^. ■ 
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B. Proof of Theorem [2] 



Assume that N < Z. Under this assumption we shall show that a 7 minimizing o (7) 
satisfies tr 7 = A^. 

Assuming the contrary, we shall find a trace class operator a > such that for 7^ = 
(1 — £:||(t||)7 + ecr and all sufficiently small e > 0, 

^(7e) < ^(7) • (65) 

The factor (1 — guarantees that < 7^ < 1 for < e < ||cr||^^. If tr7 < A^, which 

we assume, then also tr7£ < A^ for small e and fl65l) leads to a contradiction since 7 was 
assumed to be a minimizer. 

To prove fl65|) we use convexity for the homogeneous terms in the functional and 
expand the repulsion term explicitly. This leads to 

^(7.) < ^(7) + s (tr(-V=^ + l/8)a - X((t1/')) - eRi + e^R^ , (66) 

where 

(^,(r) = K(r) - I 0^^dr', 

Ri = \W\\ (^(7) + D(p^,p^)) , 

R2 = D{Pa - \\cr\\p^,p^ - \\cr\\p^) . 

Now we proceed similarly as in the proof of Proposition [T|, letting a = ai depend on a 
(large) parameter L. More precisely, we define ai by 

(x, x') = L--'/\ir/L)gir - v')x{v' / L) q-'^Xy , (67) 

with g as in (!39l) and x > a smooth function satisfying \\x\\t = 1- Asymptotically, for 
large |r|, v^^(r) ^ {Z — tr7)|r|^-'^, which is positive by our assumption. It follows similarly 
to the proof of Proposition [T] that 

tr(-V2 - V^, + l/8)aL - ) = -^^^ / |rr^X^(r) dr + o{L~') . 

It remains to show that the terms Ri and R2 are relatively small. In the proof of 
Proposition [1] and in ( H3l) we showed that WctlW = 0{L~^) and D{pfj^,pcr^) = 0{L~^), 
which implies that Ri = 0{L~^) and R2 = 0{L^^). We can then choose L large enough 
and e small enough to conclude (165|) . 

This finishes the proof of Theorem [2l 
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IV. FURTHER PROPERTIES 



A. Properties of the Minimal Energy 

Recall that E^{N) as defined in is the lowest energy of ^^^(7) under the condition 
tr7 = A^. This energy is closely related to E^{N) defined in fl35|) . 

Proposition 4. For any Z > and N > one has E^{N) = E^{N) - N/8. 

What this proposition really says is that E^{N) + is a monotone non- decreasing 
function of A^. This, in turn, follows from the fact that we can always add mass 6N far 
away from the nuclei, with an energy as close as we please to —6N/8. This was shown in 
the proof of Theorem [21 and we shall not repeat the argument. 

Proposition 5. For any Z > the energies E^{N) and E^{N) are convex functions of 
N. They are strictly convex for < N < Z . 

Proof. By Proposition |4] it suffices to consider E^{N). The convexity follows from the 
convexity of the functional. Moreover, from Theorem [2] we know that minimizers for < 
N < N' < Z have different traces, and hence different densities. The strict convexity follows 
hence from the strict convexity of D{p,p) in p. ■ 

We now prove that the energy is bounded from below uniformly in N for fixed Z. 

Proposition 6. There is a constant C > (independent of N and the charges and positions 
of the nuclei) such that for all Z > and N > 0, Ef{N) > -CZ\ 

Remark: The proof below does not use the property that 7 < 1 and this results in the 
exponent 3 , which is not optimal in the fermionic case. Without the restriction 7 < 1, the 
exponent 3 is optimal, however. 

Proof. First, let us consider the atomic case with a nucleus of charge Z located at the origin 
R = 0. We consider '?/'(x, x') = 7^''^(x, x') as a wave function in L^(]R^) and find after 
symmetrization 



Z\r 



1 



Z\r 



/I — 



1 
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By the positive definiteness of the Coulomb kernel, D{p^, p^) > 2D{p^, az) — -D(cr^, az) for 
any az- Hence 



- ^V^, - Vz{r) - Vz{r') - ^ + i 



ij) - D{(7z, crz) 



with Vz{r) = Z\r\ ^ — J |r — r'| ^az{r')dr'. We shall choose cr^ in such a way that 

- - Ivl, - Vzir) - Vzir') - + i > 0. (68) 

From this it follows that ^(7) > —D{az, crz). Actually, we shall choose az of the form 
az{r) = Z*a{Zr) for some fixed a, which yields D{az,az) = Z^D{a,a). 

To prove fIBS]) we make an orthogonal change of variables, s = (r — r')/ y^, t = (r+r')/ ^/2, 
so that the operator on the left side of fl68|) becomes 

("^^' " 7^ + i) i ("^' " '^/^O + i ("^' " " '^/^O • 

The operator in the first brackets is non- negative (see Eq. fH8|) ). Hence it suffices to choose 
a such that the operator — — AVz{{t + a)/-\/2) is non-negative of any a G M^. Note that 
Vzi't) = Z'^V{Zt) with V{r) = |r|^^ — J \r — r'\^^a{r')dr'. After scaling and translation, 
we have to prove that — — 8V{r) > 0. For this we choose a a non-negative, spherically 
symmetric function with J a dx = 1 and with support in {|r| < 1/32}. Then by Newton's 
theorem V{r) = for |r| > 1/32, and for |r| < 1/32 one has 8V{r) < l/(4|r|2), so -V^ - 
8V{r) > by Hardy's inequality (1471) . This concludes the proof in the atomic case. 

In the molecular case we proceed as follows: We recall that we are not taking account 
of the (fixed) nuclear repulsion U, and this means that we can freely place the nuclei at 
locations that minimizes the energy E^(N). We assert that the best choice of the Hj is 
one in which they are all equal and, by translation invariance, this common point can be 
the origin. The problem thus reduces to the atomic case with a nucleus whose charge is the 
total charge Z. That the optimum choice is equal follows from the fact that for any 7 
the attractive energy for nucleus j is — J p^(r)|r — Hjl^^dr and the best possible energy is 
obtained by placing all the Hj at the point R that maximizes this integral. ■ 

B. Properties of the Minimizer 

Proposition 7. Let 'y be a minimizer of (l35l) and let = {r : p^(r) > 0}. Then the null- 
space of the spin-summed density matrix, N'tT„'~^, coincides with the set of L'^i^) functions 
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that vanish identically on M^. 

Another way to say this is that if tro-7 has a zero eigenvalue then the eigenfunction 
vanishes wherever the density is non-zero. In particular, if > almost everywhere 
then is not an eigenvalue of the spin- summed density matrix trg- 7. 

Proof. Write (tro-7)(r, r') = Xlj with ipj orthonormal and < Aj < q. Then 
M.^ \ = f]j{T^ '■ = 0}, and if (p = a.e. on then obviously '~f(f = 0. Conversely, 

let V9 G ATtro- 7 and consider 

One has tr7£ = tr7 < A^, < 7^ < 1 for < £ < Ai and 

^'J' = (tr.7)'/' + v/F|^)(^|+ (v/a7^- v^) 1^1) (^1 1 • 

As noted in the introduction, it follows from convexity that minimizing £^ for density 
matrices < 7 < 1 with q spin states is equivalent to minimizing under the condition 
< 7 < 5' without spin. Hence 

E^{N) < ^(7.) = ^(7) - v^CM + Ois), 

where 

cM^jj ^'■•)-^'';>;,';><^') dr = E A // ^'■•'•*''|;'!^['''^'^'' * > . 

j 

Since 7 is a minimizer, one has C[(p] = 0, which by the positive definiteness of the Coulomb 
kernel means ipip* = a.e. for all j. Hence ip = a.e. on M^. ■ 

At the other end of the spectrum of 7, we comment on the eigenvalue 1 of the minimizer. 
Consider the minimization problem (1351) without the constraint 7 < 1, 

^boson^^^ = mf{^(7) : 7 > 0, tr7 < AT} . (69) 

This energy can be interpreted as the ground state energy of A^ bosons in the Miiller model. 
Obviously, E^°^°^{N) < E^{N) with equahty for A^ < 1. For large values of A^ we expect 
them to differ, however. 
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Proposition 8. Assume that i?<°^°°(A^) < E^{N) for some N and Z . Then any minimizer 
7 of (135!) has at least one eigenvalue 1. 

Proof. Assume, on the contrary, that 7 < 1 and let denote a minimizer for fl69l) . (The 
existence is shown in the same way as in the proof of Theorem [H) Then 7^ = (1 — 5)7 + e'jb 
satisfies tr7£ < and < 75 < 1 for sufficiently small e > 0. Moreover, by convexity, 

S^ile) < (1 - e)E^{N) + £E^°^°'^(Ar) < E^{N) , 

contradicting the fact that 7 is a minimizer. ■ 

It is not difficult to see that E^{N) ~ A^^/^Z^ for large A^ and Z, while ^^°'^°'^(A^) ~ NZ"^. 
Hence clearly ^^°'^°°(A^) < Ef{N) for large A^ and Z. 

find numerically that in fact occupation numbers that correspond 
to core electrons of large atoms all have the value one. 

Proposition 9. Let 7(x, x') be a minimizer of ^(7) for some N and let 7(r, r') = 
7(r, cr, r', cr) be the spin-summed minimizer. Then ^{r,r') is necessarily real. 

Proof. It suffices to show that 

7V2 is real. Write 7-'^/^(r,r') = A(r,r') + iB(r,r'), where A 
is real and symmetric and B is real and antisymmetric, whence iB is self adjoint. Define 
6 = A"^ — B^, noting that both A"^ and —B^ are positive (semidefinite). The kinetic and 
potential energy of 6 and 7 are equal. Moreover, the densities p-y^r) and ps{r) are equal. 
Therefore, we just have to show that the exchange terms favor 6, i.e., X(5^^'^) > X(7^/^). 

To prove this assertion use the concavity of X(-) to conclude that > X(|y4|) + 

X{\B\), where \A\ = and \B\ = ^/Wb = On the other hand X{-f^/^) = 

X(A) + X{B), with the obvious meaning that X(A) = ^ J |A(r,r')p|r — r'\~^drdr' and 
similarly for X{B). 

To conclude the proof we have to show that X{\A\) > X{A) and > X{B) if 

B ^ 0. For the first, we write A = A^ — A^ and \ A\ = A^ + where A± are both positive 
operators. Clearly, the cross term J A+(r, r')y4_(r', r)|r — r'\^^drdr' > since |r — r'| is 
positive definite. The same argument applies to iB = B^ — B_, but now we want to show 
that / 5+(r,r')5_(r',r)|r - r'|-irfrrfr' > unless B = 0. 

To show this we use the fact that the positive definiteness of the Coulomb kernel implies 
that J a(r, r')/5(r', r))|r — r'\^^drdr' is (operator) monotone in a and in /3. Therefore, it 
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suffices to show positivity for selected eigenf unctions of B±. That is, we replace -B+(r, r') 
by eigenfunctions (y9+(r)(y9+(r')* and similarly we replace i?_(r, r') by (/?_(r)(y9_(r')*. 

Since iB is imaginary and antisymmetric, however, its positive and negative spectra are 
equal, apart from sign, so B± have the same spectrum. Moreover, Bj- are complex conjugates 
of each other. Therefore, for every V5+(r) there is a (fi-ir) and the two functions are complex 
conjugates of each other. In short, it suffices to show strict positivity of / y9(r)^(y9(r')*)^|r — 
r'l^^rfr dr', but this is true as long as the function if is not identically zero (since the Coulomb 
kernel is positive definite). ■ 

Finally, we show that a minimizer of £^{"y) satisfies the variational equation fl23l) . as 
claimed in the Introduction. 

Proposition 10. Let be a minimizer 0/^(7). Then 



-|V^ - |V^, - ^,(r) - ^,(r') - -— - - 2/ij 7'/'(x,x') = ^ 2e,^,(x)^,(x')* . (70) 



Here, (p-y{r) = Vc{r) — J p^(r')|r — r'| ^dr' denotes the effective potential, /i < —1/8 is the 
chemical potential, Cj < and the ipii'^) '^'^e eigenfunctions of 'j with eigenvalue 1. 

Proof. Let /i be the slope of a tangent to the curve E^{N) at A^. Since E^{N) is convex, 
such a tangent always exists, although it may not be unique in case the derivative of E^{N) 
is discontinuous at this point. 

Since 7 is a minimizer of ^(7), its square-root 7^/^ minimizes the expression 



among all S with < S < 1, irrespective of the trace of 5^. In fact, it is even a minimizer 
if one relaxes the condition 6 > 0. This follows from the fact that X{6) < X{\6\) for any 
self-adjoint operator 6, which was shown in the proof of the previous proposition [91 

Consequently, 7^/^ is a minimizer of (I7T1) subject to the constraint — 1 < 5 < 1. From this 



J^{6) = tr (-|V' - K(r) -fx)S^ + Dips^.ps^) - X{5) 



(71) 




(72) 



""^ e=0 

The derivative can easily be calculated to be 



tr[((-iV^ 



¥^7)7 



i/2 + y/2(_iv2 



if^) - Z^- 2/i7 




(73) 
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where is defined in ( 125|) . The condition on a is that ('?/'j|c"|V'j) < for all with 
= l^j)- Hence we conclude that 

(-IV^ - v,h'/' + 7^/'(-|V2 - - - 2/i7'/' = Yl 2e#.)(^.| , (74) 

i 

with Cj < 0. ■ 

The variational equation fITOl) was obtained by varying 7^/^ instead of 7. If 7 does not 
have a zero eigenvalue (which, for a spin-invariant minimizer 7, is the case if p-y does not 
vanish on a set of positive measure, see Prop. [7]), then these variations are equivalent. Hence 
we conclude that (1701) is actually equivalent to 7 being a minimizer in case 7 has no zero 
eigenvalue. (See the discussion in Section [ID]) . 

C. Virial Theorem 

A well known property of Coulomb systems is the virial theorem, which quantifies a 
relation between the kinetic and potential energies. We state it here for an atom. 

Proposition 11. Let K = 1 (i.e., consider an atom) and let ^ he a minimizer for E^(N). 
Then 

2tr(-iV%)=tr(Z|r|-S)-I^(p„P,)+X(7^/') . (75) 

Proof. For any A > the density matrix 7a defined by 7a(x, x') = A'^7(Ar, a, Ar', a') is 
unitarily equivalent to 7 and hence satisfies < 7a < 1 and tr7A = tr7 < A^. Since 7 is a 
minimizer, the function 

^(7a) = A^ tr(-i V27) - A tr(Z|rrS) + ^ tr 7 + \D{p„ p.) - AX (7^/^) 
has a minimum at A = 1. This implies the assertion. ■ 

V. THE MULLER FUNCTIONAL AS A LOWER BOUND TO QUANTUM ME- 
CHANICS 

We are able to show that the Miiller energy E^{N) (without the addition of A^/8) is a 
lower bound to the true Schrodinger energy when N = 2, but with arbitrarily many nuclei. 
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The situation for > 2 is open. As we remark below, our N = 2 proof definitely fails when 
N >2. 

Consider the A^-particle Hamiltonian ([2]) in either the symmetric or the anti-symmetric 
A^-fold tensor product of L^(M'^,C''). For a symmetric or anti-symmetric ip we recall that 
the one-particle density matrix 7^ is defined by 

7^(x, x') = iV j V^(x, X2, . . . , XAr)^(x', X2, . . . , Xat)* ■ ■ ■ d^N • 

Proposition 12. Assume that N = 2. Then for any symmetric or anti- symmetric normal- 
ized ip, 

Proof. Since {ipl X]^=i(^i^i ~ ^c(rj))|V^) = tr(— — V^(r))7^, we have to prove that 



|ri — J 2|r — r'l J 2|ri— 

By (!5T|) it suffices to prove that for any characteristic function x of a ball (or, more generally, 
for any real-valued function x) 

2 j X(ri)|^(xi,x2)|\(r2)rfxirfx2 + j x(r)|7y^(x, x')|\(r')rfxc/x' > (^j x(r)7v,(x, x)(ix^ . 

(76) 

Introducing \1/ as the (non-self-adjoint) operator in L^(M^) with kernel \/2ip('x,x'), we can 
rewrite the previous inequality as 

tr x^^X^ + tr Xl'J\lT ^ ^^^)' • (77) 

The proof of this inequality can be found in [46i]. For completeness, we present the proof 
here. 

Note that ^^^^ = 7^, so \& = 7y^V for a partial isometry V. Since ifj is (anti-) symmetric, 
\l>tv]> = C'-y^C, where C denotes complex conjugation. Hence V^7^V = C'j^C and, since the 
square root is uniquely defined, 

VhTv = C^fc . (78) 

1 /2 

We write 5 = 7^ for simplicity and consider the quadratic form 

Q{A,C) = ^{2tTA^C6 + tTAUVCV^6 + tTVA^V^6C6) . 
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We consider this quadratic form on the real vector space of real operators, i.e., operators 
satisfying 

CAC = A . (79) 

Note that Q{A,A) = ^ (ti A^A5 + tr A^VAVU) and that, by Schwarz's inequahty, 

{tr AUVAV^6f < {tTAUA6){tiVA^VUVAVU) . 

Recalhng (178!) and (179|) we thus see that Q is positive semi-definite. This implies in particular 
that 1)^ < Q{Xj !)• This is the desired inequality fl77j) . since Q{1, 1) = tr 7^ = 2, 

Qix) = |(trx5x'^ + trx<5VxVt(5) = |(tr x<5x'5 + tr x^X^^) and 

g(l,x) = ^(Strx^' + trxVt^^V) =trx7^ • 
Here we used (1781) once more. ■ 



The obvious generalization of inequahty (17^ to > 3 is not true, as the paper [43 1 
shows. But this does not mean that the Miiller energy is not a lower bound to the true 
energy. There is some numerical evidence for this, as mentioned in A. 3 of subsection II B[ 
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